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Abstract 

The supersonic flow past a thin straight circular cylinder is investigated. The asso- 
ciated boundary layer flow (i.e. the velocity and temperature field) is computed; the 
asymptotic, far downstream solution is obtained, and compared with the full numerical 
results. 

The inviscid, linear, axisymmetric (temporal) stability of this boundary layer is 
also studied. A so called “doubly generalized” inflexion condition is derived, which is 
a condition for the existence of so called “subsonic” neutral modes. The eigenvalue 
problem (for the complex wavespeed) is computed for two freestream Mach numbers 
(2.8 and 3.8), and this reveals that curvature has a profound effect on the stability of 
the flow. The first unstable inviscid mode is seen to rapidly disappear as curvature 
is introduced, whilst the second (and generally the most important) mode suffers a 
substantially reduced amplification rate. 


1 This research was suppported by the National Aeronautics and Space Administration under NASA 
Contracts No. NAS1-18605 and NAS1-18107 while the author was in residence at the Institute for Computer 
Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665 


i 


-1 - 


1 . I nt roduct 1 on 

The current and proposed development of high speed flight vehicles 
has rekindled the general research effort Into supersonic and hypersonic 
flows. One of the key areas of aerodynamic study Is that of boundary 
layer stability/transition to turbulence. In the case of compressible 
flow, Tol lmlen-Schl Icht Ing, Cortler and Invlscld Instabilities are all 
poss I bl e . 

The problem of the stability of axl symmetric flows is of obvious 
relevance to flight vehicles, for example to the flow over fuselages 
and engine cowlings. In a recent paper Duck and Hall (1988a) used 
triple-deck theory to consider the linear (and weakly non-linear) 
viscous instability of an ax i symmetric boundary layer in a supersonic 
flow to axisymmetrlc instabilities. It was found that viscous modes 
can exist In pairs (i.e. for a given body radius, there exist two neutral 
wavenumbers with two corresponding wavespeeds), and that at a given Mach 
number such modes occur only for a body radius less than a critical 
value (dependent on Mach number). 

In a second paper, Duck and Hall (1988b) went on to consider 
non-axi symmetric disturbances. These were generally found to be more 
Important than axisymmetrlc viscous modes (possessing generally larger 
growth rates and occuring at larger body radii), whilst again it was found 

that neutral modes existed in pairs at body radii less than some 

I 

critical value (dependent on the Mach number and azimuthal wavenumber). 

However it is generally found In the case of supersonic flows that 
inviscid disturbances are more important than viscous disturbances 
(this is in contrast to incompressible flows where viscous 
instabilities are generally dominant). 

One of the earliest attempts to study inviscid compressible stability 
was made by Kuchemann (1938); in this study, the temperature gradient 
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and the curvature of velocity profile (together with the effects of 
viscosity) were both neglected, assumptions which it turns out cannot 
be properly Justified. The work that provided a key to understanding 
this type of instability was Lees and Lin (1946), in which a rational 
asymptotic approximation was developed, analogous to the incompressible 
work of Lin (1945). It was found that the quantity 

0 f 0U*1 

— x p* — x (where u* denotes the velocity tangential to the 

3y 1 8y J 

surface, y* the coordinate normal to the surface, and p* the fluid 
density) plays a key role, very similar to that of 82u*/8y*2 in 
Incompressible theory, and as such may lead to a "generalised inflexion 
point" type of instability if this quantity is zero. It was shown that 
unlike incompressible Blaslus type layers, the flat plate compressible 
boundary layer can be unstable to purely inviscid modes. This 
(two-dimensional) work on compressible boundary layer was then extended 
to three dimensions by Reshotko (1962). 

However the major differences between incompressible and 
compressible theory were not fully uncovered until extensive numerical 
calculations were possible. The first of these, by Brown (1962), was 
followed by a series of computational studies by Mack (1963, 1964, 
1965a,b, 1969, 1984, 1987). A further important difference with 
incompressible results was then revealed, namely that compressible 
theory predicts an infinite sequence of additional modes. These are 
referred to as higher modes, and are of great importance for boundary 
layers since it is the first of these (the so-called second mode) that 
is then most unstable according to inviscid theory. 

In the light of this numerical work, the prediction of Lees (1947), 
that cooling the wall acts to stabilize the boundary layer, turns out 
to be a little misleading (cooling can actually destabilise the flow, 
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according to Mack 1969, 1984, 1987); In this case, although 
the "generalised Inflexion point" of the profile may disappear wi 
cooling, these additional modes persist. 

In the light of this work on planar boundary layers, we now t 
to consider the inviscid axisymmetric stability of the boundary 1 
on a straight circular cylinder, the generators of the cylinder 1 
parallel to the flow. In particular we wish to investigate the e 
curvature plays on the stability of the flow, and so we postulate 
that (generally) the radius of the body is of the same order of 
thickness as the boundary layer. Consistent with this we choose 
prescribe planar conditions at the "leading edge" of the cylinder 
although the techniques to be described could be readily 
extended to other leading edge conditions (e.g. "rounded tips"), 
approach may be fully justified if we restrict our attention to t 
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We take the z* axis to coincide with the axis of the cylinder, 
r* the radial coordinate, and 8 the azimuthal angle. a* Is 
the radius of the cylinder, which Is taken to be Independent of both 
z* and 8. The velocity vector is taken to be v* - (v^ , Vj , v^) in 
the (r*, 8, z*) directions respectively, and T* to be the 
temperature of the fluid. Throughout we assume the flow to be completely 
independent of 8, and it Is also assumed that the azimuthal velocity 
component v^ - 0. The (full) equations (in the cylindrical system) 
of continuity, momentum and energy, are, respectively (see Thompson 1972) 


dp* 

at* 




( 2 . 1 ) 



0 p* a 2 r * r* a 2r* z * 1 

3r* 9r* 0z* r ^r r 


( 2 . 2 ) 



8p* a S 2 * r* a ^z* z* ^z* r* 
0z* + 3r* + 3z* r* 


(2.3) 


1 JL 

r* 3r* 


(K* r* |^) + (K* • (2. A) 

3r 3z 3z 


Here p* denotes the density of the fluid, p* the pressure, c p 


the specific heat at constant pressure, K* the coefficient of heat 
conductivity. The Eulerian operator is defined as 


d a * a * a 

dT* " aT* + V ! dr* + V 3 9?’ 


(2.5) 


and the viscous stress components (assuming Newtonian flow) are 


-5- 


... -* 5 ^- + j < k * - c*) ’ • **. 


Vr* ' 2 " 


(2.6) 


W ' 2 "‘ &- * 5 <>* --*>’■ Jf* 


* ★ 
3v„ 3v, 


2 r * z * ” 2 z * r * “ *** [ 9r * + 0 Z *] ' 


(2.7) 


(2.8) 


The dispersion function r* is defined to be 

r * - 2 “* K* r * * D * V * 2 + T < x * - <•*> < 7 ■ ?*>• 

where the rate-of-deformat ion tensors are 

9v l 

D i, 

r r 3r 


(2.9) 


( 2 . 10 ) 


z*z* 3z* ' 


( 2 . 11 ) 


r*z* " ? 


3v* 


3z* 

+ 

+ 3r* 


( 2 . 12 ) 


p* denotes the first coefficient of viscosity, and X* the second 
coefficient of viscosity ( - 1.5 x jt*) . 

We now go on to assume a perfect gas equation of state, namely 


P* - 


R T 


(2.13) 


(where R* is the gas constant). We also assume that n* is 
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solely a function of T* (to be prescribed later). 

The surface of the cylinder lies along r* - a*, z* > 0, along 
which we assume v* - 0. We also assume the surface of the cylinder 
is insulated, and so 


8T* 

0r* 



(2.14) 


Conditions at z - 0 must be specified. For the purposes of this 

paper, we assume that the boundary layer at z — 0 has zero thickness 

(implying planar conditions prevail); a similar assumption was made 

by Seban and Bond (1951) and their comments regarding this assumption 

are valid here. Further, since the cylinder is taken to be straight 

and thin, to leading order the far field is taken to be uniform, with 

★ 

velocity vector (0, 0, U^) . The problem is now finally closed, and In 
the following section we go on to consider the basic boundary layer flow 
on the surface of the cylinder, obtained by an approximation of the 
governing equations detailed above. 
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3 . The basic flow 

I 3.1 The boundary layer approximation 

Here we consider the steady boundary layer approximation for 
the basic flow, derived from equations (1.1) - (1.4). A fundamental 
(and important) component of this paper is the inclusion of curvature 
terms in the governing equations; we achieve this by generally taking 
the body radius to be of the same order as the boundary layer thickness 
(a similar approach was adopted by, for example, Seban and Bond 1951, 
Clauert and Lighthill 1955, Stewartson 1955, Bush 1976 and Duck and 
| Bodonyi 1986). 

With the formation of a thin boundary layer, (comparable In thickness 
to the body radius) we expect the following classical assumptions to hold 


L_» J_ , 

3r* 3z* 


and 


v 3 * » vj*, 


(3.1) 


(3.2) 


hese 

orde 

rs wi 
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D- 

(1.4) 

then 

reduce to 





3 

3r* 

(P*vj*) + ^ + 

3 

3z* 

(p*v 3 *) - 0, 



&£ 

- o, 





3r* 




■ r * 

3v 3 * 


* 3v3*1 _ 1 3 

r * 

* Bv3*] 

b 

^ - 

+ v 3 

3z* J " r* 3r* 




1 ★ 

r * 3T* A * 3t*1 

* 

3v 3 * 

i 

. i a 

* 3T* 

j c p p 

r 1 3r* + V3 3?J 

- b 

1^1 

r* 3r* 

K r — - 
3r*j 


(3.3) 


(3.4) 


(3.5) 


(3.6) 
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The last of these equations assumes Cp to be a constant. Notice that 
the pressure Is everywhere uniform on account of (3.4), together with 
our inherent assumption that the body radius is of the same order as 
the boundary layer thickness, i.e. very thin. 

One important consequence of (3.4) is that the equation of state may 
be written in the following form 


p*T* - Pco * T.*. (3.7) 

where subscript <® denotes freestream conditions. 

It is now convenient to Introduce non-dimensional quantities 
( v i . v 3 , r ,z,T ,p,fi) 


- (Re v 1 */U 0 c*,v 3 */U a) *,r*/a, Re-l Z */a, T*/T«o*. P*/Pco* . (3.8) 

where Re is the Reynolds number, based on body radius a*, namely 

ii * * * 

Re _ ^22 , (3.9) 


which must be assumed large if the assumptions (3.1) and (3.2) are to 
be valid. 

Equations (3.3) - (3.6) may then be written in the following 
non-dimensional form 

Ir (?) ♦ 3 ♦ S IP] - 

- 0. (3.11) 



(3.12) 


(3.13) 
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Here a Is the Prandt 1 number, namely 

<T - (3.14) 

which In this paper we shall assume to be a constant; y Is the ratio 
of specific heats, and the Mach number, namely 


Mo, " U»*/(7 R*Tco*)*- 


The boundary conditions are 


8t 

3? 


r-1 


- 0 , 


v l - v 3 - 0 on r - 1 , 


V3 -♦ 1 

T -* 1 


as 


r -♦ ®. 


(3.15) 


(3.16) 


To close the problem formally, we require a relationship linking 
/i to T. Here we take the simplest form, namely the linear Chapman 
law (Stewartson 1964), 


H - T, (3.17) 

although, here, conceptually, there would be no difficulty incorporating 
a more complex viscosity/temperature law. 

3.2 Numerical solution 

As z -» 0 we specify (I) the solution approaches planar conditions 
(Stewartson 1964) and so we expect the solution to become singular. 

This latter condition renders the problem in its present form 
inappropriate for numerical treatment. Instead we write 
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v i - f_1 v i (v. 0 . 

v 3 “ v 3 (V, O . 


T - T 

(3.18) 

with f - zi , 

(3.19a) 

v - (r-l)/r. 

(3.19b) 

The "hatted" functions are now expected to be 

completely regular 

as f -* 0 , approaching their planar counterparts. 

Equat ions (3 . 10) , 

(3.12), (3.13) then become, respectively 


T V 1 t) - V 1 t t, " $ >7 t v 3jj + i fTv 3r 


^ ' fTv i 

+ i t 7 t v 3 - i fv 3 T f + J + “ 0. 

(3.20) 

A A A A A 

v i v 3 tj - i v v 3 \ 3r} + i rv 3 v 3r 


- T? [''Jni + + T^V3,. 

(3.21) 

V 1 T 1 J - \ V v 3 T v + \ f Tf V 3 “ t2 (T-1)Mcc 2 ( v 3t7> 2 


T - , ' T2f 

+ a (1 V + a + a(l + r7j) T r 

(3.22) 


subject to 
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- 0 on i) - 0 , 

v l - v 3 - 0 on »j - 0 , 

T -» 1 as r) ■* «j, 

V 3 -* 1 as T) -» <*>. (3.23) 

It Is possible to define a st reamfunct I on which would ensure 
the continuity equation (3.20) Is always satisfied; however in this 
case, in addition to the order of the momentum equation being increased, 
the coefficients of the equations become considerably more complicated. 
Further, it does not appear possible to introduce a Howarth-Dorodni tsyn 
(Stewartson 1951, Moore 1951) like transformation which in the planar 
case considerably simplifies the governing equations. For this reason 
it was decided to seek a numerical solution to vj , V 3 and T directly. 
Notice setting f - 0 reduces the system (3.20) - (3.23) to the planar 
problem, namely the ordinary-differential system 

T v 1tj - V 1 T ij " iV T V 3 + i vT v v 3 ~ °. (3.24) 


V 1 v 3ij ~ \ *? v 3 v 3 tj “ j2 v 3*p? + T T t? v 3jj- (3.25) 

vi T „ - \ rjvj T, - T2 (v 3 ,)2 + T _! + T - <T,)2. (3.26) 


(again subject to (3.23)). 
The variables 


V 3 1 - v 3jJ , 


T 1 - T,. 


(3.27) 
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were Introduced, and the system (3.24) - (3.26) together with (3.27) 
were written as a system of first order ordinary differential equations, 
which were then approximated by second-order finite differences. The 
truncated system was then solved by means of Newton iteration. At each 
iteration level, the algebraic system was of block-diagonal form, with 
each block comprising 10x5 elements. 

Once the above solution was obtained, the system (3.20) - (3.23) 
was treated in much the same way, with a Crank-NIcolson approximation 
being used to approximate f derivatives (again the problem was treated 
as a system of first order equations In rf) . In this way, the solution 
was extended forwards in f. 

Fig. la shows the distribution of vj^ (t;-0) with f, and 

Fig. lb the corresponding distribution of T(tj-O). These results 
are for - 2.8, with fluid constants a - 0.72, y - 1.4. The 

(*j-0) distribution Is singular In the planar limit as f -> 0, 
and then appears to (slowly) fall continuously as f increases. The 
T(tj-O) distribution declines slightly from Its planar value at f-0. 

Results for M^ - 3.8 (same fluid constants as above) are shown 
in Fig. 2a ({"“^ 3 ( 77 - 0 ) distribution) and Fig. 2b (T(tj-O) 
distribution); these suggest the same basic characteristics as the 
lower Mach number results. 

The limit as f -» ® is of some interest. Although, as indicated 
earlier, the incompressible case in this limit is well documented, 
the particular details for the compressible case do not appear to 
have been described, although Stewartson (1964) speculates that a 
similar approach to the incompressible case is necessary. In the 
following sub-section we show this is indeed the case. 
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3,3 The far downstream flow 

We find It convenient to reconsider the system (3.10) - (3.13) 
(together with (3.17)) when studying the limit of z -♦ oo. The 
incompressible work of Clauert and Lighthlll (1955), Stewartson (1955) 
and Bush (1976) suggests that two radial lengthscales are important 
in this limit, namely r - 0(1) and r - 0(zi). Following these earlier 
works, we find It convenient to define the parameter 


e - (^ logz) -1 , 

(3.28) 

- (log 0 _1 


which is necessarily small as z -» ®. 


Guided partly by the 1 ncompress Ibl e- case , 

we expect the solution 

for r - 0(1) to take the following form 


V! - o(i/r). 

(3.29a) 

v 3 - v 3 (r, f) + 0( x /f) , 

(3.29b) 

T - f 0 (r,O + 0(l/f). 

(3.29c) 


where "barred" quantities are expected to be generally order one. 
Subsitutlon of (3.29) Into (3.10) - (3.13) Immediately reveals the 
result (neglecting 0(Vf) terms) 


a 

3F 



0, 


(3.30) 


and 


T0(7-DH» 2 c 2 



i a 

F 3F 


r 


in 

a dr 


- 0. 


(3.31) 


Integrating (3.30) yields 
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T * 

T o 3T ” ?• 


(3.32) 


where K is Independent of r. Substitution of (3.32) into (3.31) yields 


f 0 

T o 3? 




+ (7-l)K 2 H. 2 c 2 g _ Q 


(3.33) 


To facilitate the solution to (3.33) we write 


r - lnr, 


(3.34) 


and so 


f 0 ~Z 
0r 


to 2!fl 

8F 


+ K 2 £rM 00 2 ('y-l)£ 2 - 0. 


(3.35) 


This equation is further simplified by the use of a second transformation 


giving 


R - 


dr 

T 0 


8 2 T f 

8r 2 


+ K 2 f^ c 2 cr(7-1) e 2 - 0. 


(3.36) 


(3.37) 


Consequent 1 y 


T 0 - - ^2(7(7-! )€ 2 R 2 + ( A 0 R+B 0 , 


(3.38) 


where Aq and Bq are independent of R. 

Before proceding further with this solution, we consider next the 
outer solution where rj - 0(1). Guided by the above, and also again the 
corresponding incompressible results, we expect the solution to develop 


as 
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v 3 - 1 + ev 3 (ij) + 0(e 2 ) , 
T - 1 + cT 3 (i) + 0(e 2 ), 
V 1 ” + 0 ( £ 2 ) , 

where we have written, 


V 



0 


(3.39) 


(3.40) 


and also we have Implemented freest ream conditions on v 3 and T. 

if To(r,f) Is to match correctly to (3.39), we must have (to 
leading order) 

A 0 + B 0 - l* 2 Hj(y-iy<r - 1, (3.41) 

(It Is now clear that although the first term on the r ight -hand-s 1 de 
of (3.38) Is notlonally 0(e 2 ), Its Inclusion Is essential for a 
correct matching process, as Is the second term). Note that we have used 
the result 

R - r - Inr as r -» <®. (3.42) 

The matching of v 3 (r,f) with the outer solution is achieved by setting 
K - 1. 

The 0(e) corrections to (3.39) are then given by 


v 3 (»j) - 


e 


d»». 


and 


(3.43) 


Tl Cr\) - [ Aq -(T^-DM* 2 ] 


S - 


GO 


(3.44) 
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Bo t h rj and R may be regarded as forms of "optimal coordinates", as 
their use is essential for the correct matching of solutions. 

To complete the problem, we now require to specify conditions on 
r - 1. In the insulated wall case, to which the bulk of this paper is 
devoted, we require 

- 0, (3.45) 

r-1 

and so Aq - 0, (3.46) 

and this implies that 


9T 

3F 


T o| r -l - 1 + \ (7-1)Mc O 2 +0(O, (3.47) 

which impl ies 


(on account of (3.32)). 

These asymptotic results are indicated as broken lines on 
Figs. 1 and 2 for comparison with the full numerical results; the 
agreement is satisfactorily given the relative "largeness" of the 
small parameter e. 

If on the other hand the surface of the cylinder is heated or 
cooled, i.e. t|^ ^ is specified (to be T w say), then we must have 

B 0 - T w , (3.49) 

and so to leading order 

A 0 - 1 + ^(y-DM^-T*, 



r-1 1 + ?(-v-1)m2 


+ 0(c), 


(3.48) 


(3.50) 
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gl ving 



[1 + \ <t(7-1)M«, 2 -Tw] + 0(e2), 


( 3 . 51 ) 


together with, 



i- + 0(e). 

'w 


( 3 . 52 ) 


In the following section we turn our attention to the Inviscld 
instability of flows corresponding to the insulated wall class. 
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4 . Inviscid disturbances 
4.1 Disturbance equations 

We now seek to determine the effect of a small amplitude disturbance 
on the basic flow described in the previous section, to determine 
whether growth/instability can occur. We impose a disturbance whose 
wavelength is generally comparable to that of the boundary layer 
thickness and therefore also of the body radius (0(a*) ) , in which 
case the parallel flow approximation can be fully justified; this 
implies the disturbance equations are all inviscid. 

Specifically, at a fixed z station, we write 


v\* - SqO v(r) e ,Qr<Z Ct) (4.1) 

v 3 * - U ro * |w 0 (r) + «S(r)e ia(Z ' Ct) j , (4.2) 

T* - Too* | to ( r ) + «T(r) e la < z - ct >j (4.3) 

P * " P ” [W + ip(r) e i0(Z ' Ct) ], (4.4) 

p* - Poo*R*T„* 1 + «p(r)e fQ,(Z_Ct) j, (4.5) 


where i is some small (disturbance amplitude) parameter, 

t - (U„*/a*)t*, 

Z - z*/ a* , 


and so the Z - 0(1) scale is very much shorter than the z - 0(1) 
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scale, a Is the non-dimensional spatial wavenumber, and c the 
the non-dimensional wavespeed, and 

w 0( r ) " v 3 (r,z), 

T 0 (r) - T(r , z) , (4.6) 

where v 3 and T are as defined In Section 3. 

Substituting (4.1) - (4.4) into (1.1) - (1.4), (1.13), taking 
the 0(5) terms with the highest order In R yields the following 
disturbance equations 

-• + fs + T^ z + k* r + ,w ° p + 5 5 ? 4 ^ - °- ( *- 7) 


- «£ ; + + - - iL, (4.8) 

T 0 To To ^2 


ia^c 

TT 


V 


+ 


0(2 Wo V 

T 0 



(4.9) 


1_ 

T o ; 


icT + Iwq T + vTo r 


"y— 1 — — 

r + J y icp - * W 0P “ °. 


(4.10) 


together with the disturbance equation of state 


P 


T 0 P 


+ 


T_ 

T 0 ' 


(4.11) 


After some algebra the system reduces to the following two equations 
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. 2 o 

- + V w 0r v _ lp Tq-M^WQ-c)* 

r r " w 0- c L w 0- c 

i5 2 (»o-c) ; (4 .i3) 

T ° i«. 2 ' 

It was shown In the previous section that the appropriate scale as 
z -» 0 is 17 (see (3.19b)) rather than r. Consequently we find it 
convenient to use rj however, rather than r, and also to introduce 
a scaled radial velocity component ip, and wavenumber » by 

v - K<p, a - a/f 


(4.12) 


and so 


Vr\ 



*0r)<P _ *P 
"°' C 7^ 2 


Tq-H, 2 (w 0 -c) 

w 0 -c 



(4.14) 


ta 2 (wQ-c) 

T 0 


TMcc 2 


(4.15) 


Equations (4.14) and (4.15) may be combined to give 


d_ 

di? 


(wq-c)(^ + tty V ) - w 0t; * 


T 0 - Moo 2 [w 0 -c] 2 


Of 2 

T 0 


(w 0 -c)^. 


(4.16) 


This equation is very similar to the well-known planar inviscid 
equation (Lees and Lin 1946, Reshotko 1962, Mack 1984, for example), 
except for the inclusion of the single curvature term, on the 
left-hand-side of the equation; notice that allowing f -* 0 retrieves 
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the planar result . 

On rj - 0, we require 

p(Tp-O) - 0 

(the impermeability condition). The second condition is that 
is bounded as tj -> <*>. This is achieved by taking the 
r\ -» oo limit of (4.16), i.e. 

+ T+?tj ‘ (i+^) 2 " a2 I 1 - M » 2 ( 1 - c ) 2 ^- 


giving 

V ~ K 1 {±a[l-M 00 2 (1 . c) 2 ] i ( J ; + n) , , 

where K n (z) denotes the modified Bessel function of order n, 
argument z, and the sign is chosen such that the real part of 
argument is positive to ensure boundedness as tj -* <»; ^ is a 
constant. Equation (4.19) also leads to 

Vjij io <l-c)K 0 Ua[l-H. 2 (l-c)2]i(J + rj) ) 

11-^2(1-02] i 


(4.17) 


(4.18) 


(4.19) 


t he 


(4.20) 


The eigenvalue problem (for the temporal case as considered here) 
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is then, for a given a, to find c subject to boundedness as t ) -* oo, 
and such that the impermeability condition (4.17) is satisfied. 

Before carrying out a detailed numerical study of this above 
eigenvalue problem, we turn to study an important sufficiency condition 
relating to the existence of certain unstable eigensolut ions . 


5 . The doubly generalised Inflexion condition 


In the case of compressible planar flows, the existence of the 
so called generalised inflexion point, where (in our notation) 



is of great importance, as shown originally by Lees and Lin (1946) and 
confirmed by Reshotko (1962), Hack (1984, 1987) for example. 

If condition (5.1) is satsified at some point * 7 ) , then there 
exists a neutral solution, with wavespeed c, such that 

c - w 0 (»jj) , (5.2) 

provided Tg - M^ (wg-c)^ > o, (5.3) 

for all tj. 

This is a sufficient condition for the existence of a so called 
neutral subsonic disturbance, l.e. for which 

1 - < c < 1 + Vmoo" (5.4) 

using the terminology described by Mack (1984), where sonic disturbances 
have 

c - U 1 /*", (5-5) 

and supersonic disturbances have 

c < 1 " VfV (5.6) 

The condition (5.1) has a further important repurcussion, namely that 
a sufficient condition for the existence of an amplified disturbance 
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at some t) > r)\ , where r; j Is the point at which 


w 0( T i) “ 1 _ Vm*,- (5.8) 

The question then arises, what is the effect of curvature on these 
important conditions? We address this question next. We take equation 
(4.16) as our starting point, and follow the general approach adopted 
in the past to tackle inflexional instabilities arising in planar 
compressible flows (eg. Mack 1984) , although here the situation is 
more complicated because of the inclusion of curvature terms. 

Taking equation (4.16), dividing through by wq-c, and multiplying 
by ip* (where an asterisk denotes a complex conjugate) we obtain 


★ j 

a 

wq-c dr 


(w 0 -c) (if ) r + - <p r ) - W 0 ip 

r n_ 



( 5 . 9 ) 


where r - 1 + {" 17 , 


and 


X “ T 0 - K ! 0 2 (wo _ c) 2. 


(5.10) 

(5.11) 


If we now subtract the complex conjugate of (5.9) from (5.9) we 
obtain 


ip* d 

(w 0 -c)(^ r + i <p) - w 0rV p 

_ f , d 

(wq-c*) ( p * + i <p*)-w o r p* 

* 

0 

1 

0 

51 

X 

wq-c dr 

X* 


(5.12) 


After some algebra, this may be written as follows 
1 


* d? 


ip r + - <p 


* d? 


r ¥»r* + \ <P* 


- rw 


1_ <L [ w 0r l _ 1 d_ fwQ r 

w 0 -c dr [xr w 0~ c * dr [x>. 


(5.13) 
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We now focus attention on the limit of the neutral state, l.e. if 

c - c r + i c | 

then c j -» 0 . (5.14) 

We may write X* ~ X In this limit without any difficulty 
(assuming the wave is not given by either of (5.5), which we shall 
see is outside of the scope of the following, anyway). However we exert 
some care in the treatment of the right hand-side of equation (5.13), 
which we now write as 


1 d 

J V*(V’r + " >V(W* + p P*) 

2Iriv’l^ic| d 

!&£_ 

r dr 

X 

1 wq — c I ^ dr 

X r 


We now use the following arguments: ( i ) as cj -» 0, the derivative 
of the term in parentheses on the left-hand-side of (5.15) is always 
zero, except possibly at the point, rj, where wq - c; 

(ii) the term inside the parentheses must be zero at the wall (r - 1), 
and asymptote to zero at infinity if the wave is subsonic : ( i i i ) the 
r i ght -hand-s i de acts as a delta function at rj as cj -» 0, unless 


- 1^1 
dr T 0 r r-r , 


- 0 , 


or equivalently 


d_ 

dtj 



(5.16) 


(5.17) 


where r \ j - (rj-l)/f. (5.18) 

(note x( r ~ r j) ” Tq) • This condition is clearly required in order 
to avoid a finite jump in the term in parentheses on the left-hand-side 

of (5.15), and the subsequent contradiction. Equation (5.16) is a 
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"doubly generalised Inflexion condition", and Includes a curvature 
term, not present In planar studies. 

We thus see that (5.16) Is a necessary and sufficient condition for 
the existence of so-called subsonic modes. In the following section 
we carry out a numerical study of (4.16); as we shall see, (5.16) gives 
us an extremely useful guide to the behaviour, nature, and existence of 
the various modes of instability present. 
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6 . Solution of disturbance equation 

6.1 Numerical method 

For the purposes of numerical solution, equations (4.14) and (4.15) 
were chosen (In preference to (4.16)). A fairly straightforward 
Runge-Kutta scheme was applied to this system, with the shooting 
commencing at some suitably large value of rj, with (4.19) and (4.20), 
and the computation procedi ng Inwards, towards tj-0 . The Impermeability 
condition at r) - 0 was satisfied by choosing the appropriate value of 
c (by means of Newton's method). 

In a number of computations It was found advantageous to divert 
the computation below the real rj axis, In particular when iwq-ci 
was small (If lm(c) < 0 this procedure must be used). A similar 
technique has been (used by Mack (1965), using a method based on that 
of Zaat (1958). 

6.2 Doubly generalised Inflexion point results 

Before presenting details for the eigenvalue problem per se, defined 
In Section 4, we return briefly to consider further results for the basic 
flow. It was shown in Section 5 how the so called doubly generalised 
Inflexion points are likely to play an Important role in the stability 
of the flow. Consequently we return to consider the two examples 
studied In Section 2.2, namely Me - 2.8 and M* - 3.8. In particular 
we are interested In the existence of doubly generalised Inflexion 
points . 

Fig. 3a shows the axial variation of position of the doubly 
generalised inflexion points for M* - 2.8. The point f - 0 
corresponds to the leading edge of the cylinder, and as such 
corresponds to the planar case (as a result of our basic assumptions). 
There are two particularly striking features to this distribution: 

(i) that these points occur in pairs and (li) there exists a critical 
value of f, downstream of which no such points exist. The upper 
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points are an extension of the generalised inflexion point found important 
in planar cases, whilst the lower points rise off the surface of the 
cylinder tj - 0 , to ultimately merge with the upper branch at f = 0.059. 
It is remarkable how the doubly generalised inflexion points 
disappear at such a small distance downstream of the leading edge. 

It was also shown in Section 5 how neutral solutions with 
wavespeed 

c “ wo(’M) • (6.1) 

will occur, provided 

1 - Vm® < c < 1 + Vm®. (6.2) 

and so in Fig. 3b we show the axial distribution of wg(»7j) for 

“ 2.8. Because of the restriction (6.2) It Is seen that subsonic 
inflexional modes of instability will only occur for 0 < f < 0.047, 
implying that such modes will completely disappear at Just a distance 
approximately 0.0022 x Re body radii downstream of the leading edge 
(although other modes types are certainly possible); consequently in 
this case we expect this mode will disappear before the doubly 
generalised inflexion points have merged. There are certain similarities 
here with the effect of cooling of planar boundary layers (Lees 1947, 

Mack 1987), which causes a similar effect on generalised inflexion points. 

We next turn our attention to results for the higher Mach number 
considered previously, M^, - 3.8. Figure 3c shows the axial variation 
of the doubly generalised Inflexion points in this case; the general 
characteristics are the same as those of Fig. 3a, except the range of 
f for which such points exist is increased. The corresponding 
distribution of wq(t/j) is shown in Fig. 3d; this too is similar to 
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the corresponding M* - 2.8 distribution shown In Fig. 3b. In the case 
of Moo “ 3.8, Fig. 3d Indicates that neutral subsonic Inflexional inodes 
will disappear a distance approximately 0.013xRe body radii downstream of 
the leading edge. 

Guided by the above observations, we now turn our attention to the 
eigenvalue problem for the two cases Moo — 2.8 and Moo — 3.8. 

6.3 Growth rate results 

Figure 4a shows the variation of cj with a (where c - c r + 1 cj), 
for the case Moo- 2 . 8 , at f -0 (and hence corresponds to a planar 
example). The corresponding results for c r are shown In Fig. 4b. Here, 
and in all results to follow, neutral points are denoted by a cross. 

These results (which are typical of previous planar results - see for 
example Mack 1987) show two distinct unstable modes. The first (to be 
referred to as mode I) originates as a sonic neutral disturbance (with 
cj - 0 , c r - 1 - Vh®) at cr - 0 , and terminates as a neutral Inflexional 
subsonic mode at a a 0 . 1 , where c r — w(ij - 17 j ) a 0 . 66 ; this mode in 
fact continues, becoming a decaying mode, with cj < 0 , although we 
shall mainly concentrate our attention on growing/neutral modes). 

The second mode (to be referred to as mode II) originates at 
a a 0.4 as a subsonic neutral mode with cj — 0, c r — 1 (this may 
be regarded as a special case of an inflexional mode, with the 
generalised inflexion point occuring in the freestream). This 
mode then terminates at a a 1.13 as a (second) neutral subsonic 
inflexional instability (and at values of a greater than this 
value continues as a decaying mode, with cj < 0 ) . 

Although other modes of instability undoubtedly exist at this Mach 
number, these have considerably smaller growth rates than modes I and 
II shown here, and are subsequently much less important from a practical 
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point of view. Note that since the (temporal) growth rate Is arc], 
mode II Is the most Important. 

We now turn to results Incorporating the effects of curvature. 

Fig. 4c and Fig. 4d show distributions of c r and cj (respectively), 
with a (for M® - 2.8), at f - 0.02. Although the qualitative 
features resemble those of the t - 0 case, the maximum amplitude of 
the growth rates is seen to be considerably reduced (in spite of the 
smallness of f) , particularly that of mode 1. 

Moving further down the axis of the cylinder, to f - 0.04, Figs. 

4e (c| distribution) and 4f (c r distribution) Indicate that mode I has 
practically disappeared, whilst the maximum amplitude of mode II Is now 
significantly diminished, terminating (at a subsonic Inflexional neutral 
point) at quite a large value of or (* 2.65), although over much of 
the range of a, this mode has exceedingly small growth rates. 

Following our comments In the previous subsection regarding wo(i 7 j) 
dropping below 1 -Vm°°> we *xpect mode I to completely disappear at 
f a 0.047 for this choice of M^,. As a consequence, the next set of 
results (at f - 0.05) presented In Fig. 4g (cj) and Fig. 4h (c r ) shows 
just mode II, which Itself exhibits a further reduction In growth rate. 
This mode still originates as a neutral mode with c r - 1 (at a =.45); 
unfortunately our computations did not indicate a clear neutral 
solution at an upper value of a. This was due to the exceedingly small 
growth rates encountered, which were typically 0(10“^®), and hence 
were comparable with the round-off associated with the computation. (In 
the regime of larger a and very small growth rates, It was 
found to be most essential to deform the Integration contour in the 
numerical scheme, as described In Section (6.1), in order to maintain 
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numerlcal accuracy.) If a neutral point exists, as seems likely, it must 
be of the neutral supersonic class (cj - 0, c r < 1 - 1 /mco) because of 
the absence of any doubly generalised inflexion points at this value 
of f. 

As a final example of the Me - 2.8 flow, we show results for 
? - 0.2 in Figs. 4i (cj) and 4j (c r ) . These indicate qualitative 
similarity with the previous set of results; however the maximum growth 
rate is reduced by approximately an order of magnitude. Again, 
unfortunately, positive identification of an upper neutral point was 
not possible, due to the difficulties with tiny growth rates 
encountered at larger values of a. We conclude, however , that 
curvature has Important (and profound) effects: (i) annihilation of mode I 
and (ii) substantial reduction of the growth rate of mode II (although 
the range of a over which this unstable mode exists is increased 
quite significantly). 

We next turn our attention to results for Me - 3.8, and Figs. 5a 
and 5b show cj and c r distributions (respectively) with a, for 
the particular case f-0. This corresponds to the planar case as 
computed previously (Mack 198.7 for example) and thus provides a 
useful check on the accuracy of the present overall scheme, (which 
is seen to be entirely satisfactory). When compared with the 
corresponding M*, - 2.8 results (Figs. 4a, 4b), the importance of 
mode II is seen to be significantly increased (although the growth 
rate of mode I Is increased also). Generally, the ^-3.8 
distributions qualitively resemble the corresponding M^ - 2.8 
distribut ions. 

At f - 0.05 (with M* - 3.8), we see in Figs. 5c (cj distribution) 
and 5d (c r distribution) there is an approximate halving of the 
maximum growth rate, when compared with the f-0 results. Further 
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downstream, at f - 0.10 (Figs. 5e , 5f) mode I has almost 
disappeared, whilst mode II has suffered a further depletion of 
maximum growth rate. 

From our observations made in Section (6.2), we expect that mode I 
will completely disappear at f » 0.013 (where c r - wq(ij - »jj) - 
1 ■ Vmot) I this will also be an important location for the upper 
neutral point of mode II, which for f < 0.013 is of the subsonic 
inflexional kind. As noted earlier, for the previous M* considered, 
growth rates in this regime were extremely small, and so no firm 
conclusions on the behaviour of this mode in this region were possible. 
Fortunately, although the growth rates in this critical region at 
Hx, - 3.8 are small, they are nonetheless significantly larger than 
at the lower Mach number. 

Figure 5g shows the local variation of cj with a at f - 0.112 
(just below the critical value). Mode II is clearly seen to become a 
damped mode at as 0.85, with I c | l reaching a maximum at as 0.93, 
and then decreases. Unfortunately, no firm conclusions are 
possible regarding the ultimate behaviour of cj at large 
values of a, due to the smallness of I c j | . 

Figure 5h shows the distribution of C| in the same critical 
region of a, at f - 0.114 (slightly above the critical value of 
f) . It now appears that the growing mode II now terminates at 
a s 0.94, as a supersonic neutral mode, and does not continue as a 
damped mode . Instead, a further (supersonic) neutral mode has already 
appeared (at a = 0.88) and this is then the origin of a damped mode, 
which has a maximum value of |cj) at a = 0.95; Icj | then decreases, 
towards zero, and again because of its ultimate smallness, no 
conclusions regarding its behaviour at larger values of a are possible. 
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Addit ional ly, there was some evidence of another unstable mode, 
beginning ai'ound a * 1.52, but because of Its very small growth rate 
It Is not possible to be completely categorical about this; its 
growth rate was also too small to be seen on the scale of Fig. 5h. 

Thus In this case we see the presence of possibly three (supersonic) 
neutral points In this regime. 

Figure 51 details the variation of cj in the crucial a region, 
for the location f - 0.116. Again, (growing) mode II is seen to 
terminate at a supersonic neutral point, this time located at a w 1 . 

A further supersonic neutral point exists, originating at aw 0.9 
which then provides the start for a decaying mode; when compared to the 
corresponding f — 0.114 mode, the decay rate of this particular mode 
Is seen to be reduced. Further, this mode now seems to terminate 
at another neutral point (at a w 1.48). Yet another neutral point 
exists at as 1.40, which then provides the start of a second 
unstable mode (although the growth rate of this mode was so small to 
be bearly visible on Fig. 51); a total of four supersonic neutral 
points are thus observed In this region of a. 

Figure 5j shows the cj distribution at f - 0.118 In the same 
general region of a. When compared with the previous results, a 
further change to the qualitative picture is seen. Here, the original 
mode II has merged with the second unstable mode. Just two neutral 
points remain In this region, at aw 0.93 and a w 1.35, which are 
associated with the start and the terminus of the decaying mode (which 
generally has a significantly reduced icji compared to the 
previous results). The ultimate behaviour of the growing mode with 
a remains unclear, due to the smallness of icjl. 

It is interesting to note that when two modes were present (for a 
given value of a), both modes had values of c r that were 
practically Indistinguishable. Further, there was no difficulty in 
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calculatlng accurate values of c r , even at large values of cr. 

Moving further downstream to f - 0.2, the decaying mode has 
disappeared completely, and cj and c r distributions over a broad 
range of a are shown in Figs. 5k and 5C respectively. When compared 
with the Fig. 5e results, the mode II growth rates are quite 
appreciably reduced; the ultimate behaviour of lcj| at large values 
of a remains unclear. 

Further downstream st i 1 1 , at f - 1.0, results (cj shown in 
Fig. 5m, c r shown In Fig. 5n) are qualitively similar to those at 
f - 0.2, except that the maximum value of cj is significantly 
diminished, and occurs at a rather larger value of cr (as does the 
origin of this mode which occurs at cr = 0.72, compared with a = 0.4 
in the case of f — 0.2). The larger a behaviour of this mode 
is again unclear, due to the reasons described above. 

In the following section we turn to consider a number of general 
conclusions and points raised by this work. 
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1, Conclusions 

In this paper the supersonic flow over a thin straight circular 
cylinder has been Investigated. The basic boundary layer flow has 
been obtained, and the Inviscid stability of the flow has been 
studied. A condition on the basic flow for the existence of so called 
subsonic Inflexional neutral inodes of instability has been derived, 
and is found to be an extension of the generalised inflexional 
condition relevant to planar flows. 

The effect of body surface curvature is seen to immediately (and 
significantly) reduce the Importance of the "first mode" of inviscid 
Instability, which is seen to completely disappear at what could be 
a comparatively short distance down the axis of the cylinder (by 
about 0.0022 x Re body radii at - 2.8, and by about 0.013xRe 
body radii at M^ - 3.8). 

The maximum growth rate of the "second mode” of Inviscid 
instability also suffers a substantial reduction at locations 
increasingly further down the axis of the cylinder, although the 
evidence is that it does not disappear completely. 

There are certain similarities here with the effect of cooling 
a planar boundary layer (Mack 1984, 1987, for exmaple), which can also 
cause the first mode to disappear completely (cooling also causes the 
formation of a second generalised inflexion point, which with a 
progressive reduction in wall temperature eventually coalesces 
with the first generalised inflexion point). However the effect of 
cooling is to increase the amplification rate of the second mode (in 
contrast to our results featuring curvature). 

It is particularly interesting to note that although inviscid 
modes of instability are generally regarded as more important/ 
unstable than viscous modes of instability in the case of supersonic 
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planar flows., the work of Duck and Hall (1988a,b) on viscous 
axisymmetrlc flows indicates that a reduction In body radius 
(equivalent to a further downstream location In our context) 
causes an Increase In amplification rate. Thus It is entirely 
possible that with axisymmetrlc flows, regimes may exist where viscous 
instability is dominant. 

It Is to be hoped that this study will provide a basis for the 
study of flows over further and more practical geometries, such as 
cones. One Important omission to the physics of that problem, which 
must ultimately be resolved, Is the exclusion of any shock waves in 
the basic flow (however this may be justified by the restriction of 
thinness), although previous works on planar flows (cited throughout 
this paper) all have this same omission. It is also 
hoped that these results will provide material for comparison with 
finite Reynolds number computations. 
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Fig 2a Variation of £ v 3tJ (r>-0) with 
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Fig 2b Variation of T(ij-O) with 



Fig 3 a Variation of transverse positions of inflexion points ( 17 j ) 
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Fig 3b Variation of VQ(ri-r}[) with 
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Fig 3d Variation of 'MQ(rf-rji) with 
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Fig 4a Variation of cj with a, M^, - 2.8, f - 0 (pi anar case) 



Fig 4b Variation of c r with 







Fig 4d Variation of c r with a, M*, - 2.8, K - 0.02 
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Fig 4f Variation of c r with a, Moo “ 2.8, f - 0.04 
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Fig 4g Variation of c { with a, M*, - 2.8, f - 0.05 




Fig 4i Variat ion 
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Fig 5a Variation of c 5 with a, M* - 3.8, f - 0 (planar case) 
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Fig 5c Variation of c { with a, H*, - 3.8, f - 0.05 
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Fig 5e Variation of c< with 
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Fig 5h Variation of cj with a, M* - 3.8, f - 0.114 
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Fig 5k Variation of cj with 




Fig 5C Variation of c r with 
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Fig 5n Variation of c r with 
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